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Abstract 



In this study, the internal energy and the specific heat of the one- 
dimensional Ising model obtained in the frame of the generalized statistical 
thermodynamics (GST) by R.F.S. Andrade (ref.[19]), have been extended 
somewhat and compared with the internal energy and the specific heat which 
have been calculated in the standard (conventional) statistical thermodynam- 
ics. In the q ^ 1 case, approaching of the generalized internal energy and 
the specific heat to the standard well-known results, have been illustrated 
by numerical applications. In addition to these, A >> 1 case has been also 
investigated. 

I. INTRODUCTION 

Boltzmann-Gibbs (BG) statistics is essentially derived from the standard Shannon en- 
tropy 
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which is an extensive and concave quantity (the sub index 1 of 5* will be defined afterwards). 
Although BG statistics is suitable for handling a large number of physical systems, it fails 
whenever 

• the range of the microscopic interactions is large compared to the hnear size of the 
macroscopic systems (long-range interactions) 

• the time range of the microscopic memory is large compared to the observation time 
(non-Marcovian processes). 

These kinds of violations are met for a long time in gravitation [1] , magnetic systems [2] , 
anomalous diffusion [3] and surface tension problems [4] . The way out from these problems 
seems to be Nonextensive Statistical Thermodynamics which must be a generalization of 
the BG statistics in a manner that allows a correct description of the nonextensive physical 
systems as well. The study of nonextensive formalisms is very vivid nowadays in Physics and 
these formalisms keep growing along two apparently different lines: Generalized Statistical 
Thermodynamics (GST) and Quantum Group- like Approaches (QGA). Although they seem 
to be very different from each other, recently it is pointed out that there exist a connection 
between them [5]. Here, the former one will be the subject of interest and the latter is 
completely out of the scope of the present study. 

GST has been proposed by C. Tsallis in 1988 [6]. The multifractal inspired entropy 
which has been expressed by Tsallis in the form 



is a generalized entropy having the properties of positivity, equiprobability, concavity and 
irreversibility. Furthermore, the standard additivity property has also been appropriately 
generalized (in the g 7^ 1 case it becomes non-extensive). In eq.(l), /c is a positive constant, 
q (Tsallis g-index) is a real constant associated with the number of measurement in the 
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multifractal theory and W is the total configuration number. Commencing with the entropy 
definition of eq.(l), a canonical distribution 

^9 i=l 

could be obtained for a system in a heat bath, by maximizing the entropy under suitable 
constraints [6]. 

Soon after Tsallis' work, the generahzation of the conventional concepts in various fields 
of Physics have rapidly proceeded. Among the successfully tackled concepts in the frame of 
this generalization, generalized thermostatistics [7], Boltzmann i7-theorem [8,9], Ehrenfest 
theorem [10], the distribution function of the classical and quantum gases [11,12], Langevin 
and Focker-Planck equations [13], Bogolyubov inequahty [14], fluctuation-dissipation theo- 
rem [15], von Neumann equation [16], classical equipartition principle [17], Callen identity 
[18], Ising chain [19,20], paramagnetic systems [21], infinite-range spin- 1/2 Ising ferromagnet 
[22], mean- field Ising model [23], the generalized transmissivity for spin-1/2 Ising ferromag- 
net [24], self-dual planar Ising ferromagnet [25] and Planck radiation law [26,27] could be 
enumerated. 

On the other hand, this generalization has also been succsessfuUy used to overcome the 
failure of BG statistics in some physical applications such as stellar polytropes [28], the 
specific heat of the unionized hydrogen atom [29] , Levy-hke anomalous diffusion [30] , d — 2 
Euler turbulence [31], solar neutrino problem [32] and velocity distribution of galaxy clusters 
[33]. 

Up to now, a very large number of works, which have been performed in the frame of 
GST, are devoted to the investigation of the magnetic systems. 

As an example of these, the distribution function of the paramagnetic systems might be 
noted. In this study, for the paramagnetic systems exhibiting multifractal structures, the 
connection between the Tsallis g-index and the fractal dimension and the range of the values 
of the g-index have been established [21]. 

In another work [23], using the generalized Bogolyubov inequality, the free energy and 
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the magnetization of mean-field Ising model have been obtained within GST and it is noticed 
that the range of the q values are consistent with the one obtained in ref.[21]. Moreover, 
the expression obtained for the critical temperature is exactly the same as the one found in 
ref.[18]. 

Another example for the investigation of the magnetic systems in the frame of GST is 
the Ising chain which is the subject of the present work. In his first study relevant to this 
subject, Andrade calculated [19] the internal energy and the specific heat of the Ising chain 
by making use of the distribution given in eq.(2) which had been given before by Tsallis. 
Later on, in his second investigation [20] he reiterated his calculations using the distribution 

. (3) 



where the partition function, Zq is 

w 

Z,^Y.[^-(5{l-q)s,]T^^ . (4) 
1=1 

At first sight, this distribution seems to be very similar to eq.(2), only with a small redefini- 
tion of g, however, this is not the case since this change originates from the fact that in the 
definition of the second constraint used for maximizing the Tsallis entropy, g-expectation 
values have been used. 

Here, our motivation is to improve Andrade's first work [19] somewhat by considering 
r = 3 case and to strengten Andrade's paper by comparing the obtained results with the 
standard {q — 1) ones. In the present study, the generalized internal energy and specific heat 
{q 7^ 1) have been compared with the standard results. In the q ^ 1 case, the transformation 
of the results obtained by GST to the standard results, have been observed. Moreover, by 
taking into account the N » 1 case, the approaching of the internal energy expressions to 
a constant value and thus vanishing of the specific heat has also been illustrated. In this 
manner, the present study supports the investigations of Andrade as well as emphasizes the 
fact that GST is one of the most favourable tools for the examination of the multifractal 
structures. (It is obvious that, the other studt of Andrade [20] could also be tackled in the 
same manner.) 



II. ONE-DIMENSIONAL ISING MODEL 



In the model established by E. Ising for the first time with the objective of explaining 
ferromagnetism [34], a two- valued spin variable Si taking the values of +1 or —1 has been 
allocated for each lattice site. The model Hamiltonian in one-dimension, in the absence of 
the magnetic field, could be written in the following form 

N 

- jY. s,Si+i (5) 

i=l 

where J is the exchange parameter and N is the number of total lattice sites. 

The analytical solution of the model with the transfer matrix method of standard statisti- 
cal thermodynamics (SST) using the Hamiltonian of eq.(5) leads to the following expressions 
for the internal energy and the specific heat respectively [35] : 

^ = -iVtanh(/5J) (6) 

and 

^^Np^J^sech^pj) . (7) 
Kb 

III. OBTAINING SOME OF THE THERMODYNAMIC QUANTITIES OF THE 

ISING CHAIN USING GST 

The most prominent difficulty met in the solution of systems investigated with GST is the 
determination of the partition function. This situation originates from the non-extensivity 
property which arises for the case q ^ 1- Andrade wrote the Hamiltonian of eq.(5) with a 
parameter E which calibrates the ground state energy of the system for the Ising chain in 
the following form 

N 

n = E-jY,SiSi+i (8) 

i=l 

and by taking E = —NJ such that the values of the energy levels would become nonpositive 
from —2NJ to 0, he expressed the partition function of the system in a solvable form within 



the frame of GST by combinatorial method. As a consequence, he estabhshed the following 
expression [19] : 

Er Y^i — m ( §J_\^ * H \^m+l a i\fr+u—m—s 
^ _ _ J m=Ql^s=Q \ r) s\m.\{r-m-s)\ ^u=Q Sm+l,u^'' 

-r-^r ^T-m f §J.Y * ^- V"* AJr+u-m-s 

Z^m=0 Z^s=0 \ r ) s\m\{r-m-sy. ^u=0 Sm,u-^^ 

where r — l/{q — 1). 

By writing down the definitions 

u 

Cmu = Yl Cm{n)-fm-n,u-n{-'i-T (10) 
n=0 



Cm(n) - • (11) 

n!(m — ny. 



lmn^J2Vns (m,n,s>0) (12) 

s=0 



Vns"' = +{s-m- 1X,_, + (13) 

the internal energy values for r = l(g = 2), r = 2{q = 1.5), r = 3{q = 1.33) from eq.(9) could 
be written down as : 

BJ^N , , 



iiJ'^N (\i5JN + l) 



= -7VJ - 



J7V [/? J(Ar + 1)] + 1 
/3 J27V (7V2 + TV - I) + 1/3 JiV + 1 



\i5JN 



+ 1 



(16) 



(7v+f)+/3j(7v+i); 

The variations of the internal energy with respect to the reduced temperature (1//3J) for 
eqs.(14)-(16) of GST and eq.(6) of SST have been illustrated in Fig.l. In accordance with 
the expectations, for increasing r values (g values decreasing to 1) the behaviour of the 
internal energy curve approaches to that of the analytical solution of SST. Furthermore, in 
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N » 1 case, as it is illustrated in Fig. 2, Er approaches to a constant value which is an 
unexpected result. 

From the internal energy expressions given by eqs.(14), (15) and (16), taking r = l{q = 
2),r — 2(q — 1.5), r = 3{q — 1.33), the following conclusions could respectively be obtained 
for the specific heat : 

k {l + l3JNf ^ ' 

C2 _ [4 + AfiJN + p'^.PN{N - 1)] 

T ~ [4 + AI3JN + (3^jm{N + l)f 

C3 _ 9/?V^ (A^^ + - 2N) - {N^ + 3N^) - 18/?V^ + 108p JN + 81] 

T ~ [/?3J3 {m + + 9/?2 J2 {m + N) + 27] 

The variation of the specific heat with respect to the reduced temperature for = 1 has 
been illustrated in Fig. 3. It has been noticed that, the curves plotted using eqs.(17),(18) and 
(19) of GST for r = 1,2 and 3 approaches to the curve given by eq.(7) for increasing r values 
(while q values decreasing to 1), which has been established by the transfer matrix method 
of SST. Furthermore, in >> 1 case, approaching of Er to a constant value naturally 
causes to vanish. 



(18) 



■ (19) 



IV. CONCLUSIONS 

In this study, the expressions that have been obtained in GST for the internal energy 
and the specific heat of the Ising chain have been compared with the corresponding relations 
in SST. q could be regarded as a parameter belonging to the internal energy and the specific 
heat. It has been clearly shown by numerical applications that: the internal energy and the 
specific heat established for the Ising chain in the frame of GST approach to the conventional 
results obtained by the Boltzmann-Gibbs distribution in the g — > 1 (increasing r) case. This 
conclusion conforms our expectations such as the transformation of the Tsallis entropy 
(eq.(l)) to the Shannon entropy. The parametric effect of q has been brought to a concrete 
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basis with the aid of the plots. The verification of the internal energy and the specific heat 
expressions brought in by Andrade have been supported by the standard results obtained by 
the transfer matrix method. In addition to these, the A'" >> 1 case has also been considered; 
the explicit approaching of the internal energy to a constant value which inturn leads to the 
vanishing of the specific heat has been graphically illustrated. 
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FIGURE CAPTIONS 



Figure 1. The variation of the internal energy with the reduced temperature ior N — 1. 
Figure 2. The variation of the internal energy with the reduced temperature for >> 1, 
namely N — 30. 

Figure 3. The variation of the specific heat with the reduced temperature for = 1. 
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Figure 2 
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Figure 3 
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